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Fundamental approaches to a mathematical description of the liquidus surface in silicate systems are pro¬ 
posed, which can be used in developing software for operational control of the process parameters in silicate 
production. 


Glass and glass ceramic materials based on basalt and 
similar rocks find increasingly wide application. In particu¬ 
lar, basalt fiber is distinguished by high strength, high chemi¬ 
cal resistance, and other valuable properties, which makes 
this fiber indispensable in several sectors of industry. It 
would be convenient to use industrial mining waste for this 
purpose; however, such waste is characterized by a high he¬ 
terogeneity of its composition and properties. We have de¬ 
monstrated that the following levels of natural heterogeneity 
in this waste can be distinguished depending on the type of 
effect on the technological processes: 

- the crystalline-chemical level determined by differ¬ 
ences in the chemical composition and properties within a 
crystalline individual; such heterogeneity can influence the 
concentration process but has no effect on the production 
technology of construction materials; 

- the mineralogical level related to different properties 
within the same rock; the negative effect of such inhomoge¬ 
neity can be eliminated by preliminary preparation of materi¬ 
als (homogenization); 

- the petrographic level determined by variations in the 
properties of a mineral in different parts of the bed; to de¬ 
crease the negative effect of such heterogeneity, special tech¬ 
nical measures such as averaging are required; 

- the geological level determined by the presence of dif¬ 
ferent grades and types of a mineral within the same deposit 
or ore field; elimination of the negative effect of such hetero¬ 
geneity calls for the development of special technologies. 
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The industrial mining wastes studied by us had heteroge¬ 
neities at all levels, and whereas the effect of heterogeneities 
at the geological and petrographic levels can be eliminated in 
the course of primary preparation of materials and control¬ 
ling the material quality, the negative effect of inhomo¬ 
geneity at the mineralogical and especially at the crystal¬ 
line-chemical levels cannot be fully eliminated. In this con¬ 
text, operational control of the technological process is re¬ 
quired. 

The temperature of complete batch melting is of special 
significance in the production of glass and glass ceramics. If 
unmelted crystalline phases are contained in incompletely 
melted glass, local stresses arise in the finished product. The 
presence of even a small quantity of crystalline phases may 
significantly deteriorate the external appearance of finished 
articles and decrease their mechanical strength. 

In using traditional materials, the batch composition and 
the process temperature are selected in such a way as to 
avoid the emergence of the crystalline phase in the course of 
glass melting. However, this condition cannot always be sat¬ 
isfied in producing glass or glass ceramics from mining 
waste due to the unstable composition of the material recy¬ 
cled. For the operational control of the process, it is neces¬ 
sary to assess the degree of the effect of variations in the 
mineral and chemical compositions of the material on the 
temperature of complete melting. 

Most binary and many ternary systems within the limits 
of the basaltoid melt field (CaO - MgO - A1 2 0 3 - Fe 2 0 3 - 
Na 2 0 (K,0) - Si0 2 system) are well studied and represented 
as phase diagrams [1,2]. However, this form of representing 
the liquidus temperature is not very convenient for current 
practical control of the batch preparation process and the 

80 

0361-7610/02/0304-0080$27.00 © 2002 Plenum Publishing Corporation 





Programming and Calculation of Process Parameters in Production of Silicate Materials 


81 


melting and working temperature. With respect to binary and 
ternary systems, this form allows for fast determination of 
the liquidus temperature for any glass composition, but with 
significant errors. However, in the case of a four-component 
system, such determinations are neither accurate nor fast. 

Methods for mathematical description of property-ver- 
sus-composition dependences for multicomponent systems 
have been developed in recent years [3 - 6]. This problem 
can be relatively easily solved for simple binary eutectic sys¬ 
tems. In particular, I. B. Bobylev and V. N. Anfilogov [7 - 9] 
proposed a method for the calculation of the liquidus line for 
oxide - orthosilicate binary systems using the Schroeder - 
Van Laar equation: 


lnJV = 




R, 


where N is the molar part of the component crystallizing in 
the particular segment of the phase diagram; A H m and T m 
are, respectively, the enthalpy and melting temperature of the 
particular component. 

According to the data of these authors, silicate melts, in 
accordance with a polyelectrolyte model, are semi-ideal mix¬ 
tures, and the above equation can be applied to them without 
replacing the concentration with the activity. However, such 
calculation involves certain difficulties even for the simplest 
cases. This is related to the fact that melting of orthosilicate 
is accompanied by its partial decomposition into a metal 
oxide and more complex polymer forms. Furthermore, data 
on melting enthalpy are missing for many compounds, even 
for binary systems. Therefore, the authors observe that the 
main advantage of the proposed calculation scheme is not 
knowledge of the liquidus temperature but the possibility to 
determine the concentrations of simplest polymeric forms in 
silicate melts. 

Vast opportunities are provided by a mathematical de¬ 
scription of the liquidus line or the liquidus surface for bi¬ 
nary, ternary, and multicomponent systems based on experi¬ 
mental data [3], Such calculations for simple eutectic sys¬ 
tems do nor present significant difficulties. Ternary systems, 
for instance, X-Y-Z systems, with a small exception are 
approximated by equations of the following form: 

T z — Ax 2 + By 2 + Cxy + Dx + Ey 


or 


T z = Ax 2 + By 2 + Cxy + Fx 2 y+ Gxy 2 + Dx + Ey, 

where T z is the estimated temperature within the bounds of 
the crystallization field of the phase z; x and y are the weight 
contents of the phases x andy, %. 

For instance, the system CaO - MgO - A1 2 0, - FcyO, - 
FeO - Na 2 0 (K 0 O) - Si0 2 is characterized by a large quan¬ 
tity of intermediate two-, three-, and four-component mineral 
phases and binary and triple eutectics; consequently, the 


problem of a mathematical description of the liquidus sur¬ 
face within the limits of this system is difficult to solve. At 
the same time, a basaltoid composition has a relatively nar¬ 
row variation interval of the Si0 2 content; therefore, in most 
cases basaltoids do not go beyond the bounds of the 
pyroxene - feldspar (quartz) - (magnetite) field or the albite 
NaAlSi 3 0 8 - anorthite CaAl 2 Si 2 O g - sodium silicate 
Na 2 Si 2 0 5 - wollastonite CaSiO, - ferrosilite FeSi0 3 - ensta- 
tite MgSi0 3 - silica minerals (quartz, tridimite, cristobalite) 
Si0 2 field. Therefore, we restricted ourselves to describing 
the liquidus surface within a sequence of the most abundant 
ternary (and pseudotemary) systems that do not go beyond 
the bounds of the specified field. 

The successive approximation method was used for the 
mathematical description of the system. The description of 
the liquidus surface in a ternary system A— B - C started 
with a description of the liquidus lines of the binary systems 
A - B, A - C, and B - C. The liquidus line for the A - B bi¬ 
nary system on the segment from A to the double eutectic 
was described by the equation 

t=t a +k b\ B + k b2 Bl+ ••• +k bn B "’ (!) 

where t is the estimated liquidus temperature, °C; t a is the 
melting point of the mineral phase A, °C; B is the weight 
content of the phase B, %; k hl , k b2 ,k hn are the 
dimensionless coefficients. 

The liquidus line on the segment from B to the eutectic 
was described by the equation 

t = t b + k al A + k a2 A2+ ■■■ + Kn A '^ 

where t b is the melting temperature of the mineral phase B, 
°C; A is the weight content of the phase A, %. 

In the case of a simple eutectic system, a higher tempera¬ 
ture value is true for any point of the composition. It was first 
assumed that each branch of the liquidus line can be de¬ 
scribed by a first-order equation. The regression equation co¬ 
efficients were calculated using the least squares method. If 
the maximum deviation of the estimated temperature values 
from the experimental values did not exceed 10°C, it was 
considered that the equation adequately describes the 
liquidus line. 

With larger deviations, the liquids line was described by 
a second-order equation. Then the procedure was repeated. If 
the discrepancies between the experimental and the esti¬ 
mated data again exceeded 10°C, the liquidus line was de¬ 
scribed by a third-order equation, etc. In most cases the sec¬ 
ond-order equations yielded satisfactory convergence. 

In rare cases where the deviation of the estimated data 
from the experimental values in the description of the 
liquidus line was not of a specific pattern, trigonometry func¬ 
tions were used. If the crystallization field of a third phase 
was observed within the limits of a binary system, the 
liquidus surface within this field was described by an addi¬ 
tional equation. 
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After the description of the liquidus lines in the binary 
systems is completed, it is possible to pass to describing the 
liquidus surface in ternary systems. The equation describing 
the liquidus surface for the A - B - C system within the li¬ 
mits of the crystallization field of the component A will have 
the following form: 

t=t a + k bl B + k b2 B2+ ■■■ +k bn B " + k c\ C + 

k c2 C 2 + ... + k cn C n + k bcl BC + k hc2 B 2 C + 

k bd BC 2 + ... + k hcm B n C n . (2) 

Comparing Eqs. (1) and (2), one can conclude that 

At = k bcX BC + k bc2 B 2 C + k hc2 BC 2 + ... + k bcm B » C'\ 


mial. The lower boundary of the confidence probability is 
equal to 0.9 (or 90%), and the minimum number of degrees 
of freedom is equal to 4. Assuming a degree of the polyno¬ 
mial equal to 1 (i.e., a linear dependence), the minimum 
number of experimental points is equal to 6. As the degree of 
the polynomial increases, the number of degrees of freedom 
decreases, which increases the confidence interval. 

Thus, with confidence probability 0.9 and the number of 
the degrees of freedom equal to 4, the confidence interval 
will be 

b: ± 1.066 — , 


where At is the difference between the liquidus temperature 
(or the experimental temperature) in a particular point not 
lying on the side of the triangle ABC and the estimated value 
of this parameter calculated from the binary system equa¬ 
tions. 

Having calculated the latter value for each composition, 
for which the experimental temperature of complete melting 
is known and having calculated At, it is possible to determine 
the coefficients k bcl , k bc2 , k bc2 ,..., k bcm . Similarly to the pre¬ 
vious cases, this calculation was made by the successive ap¬ 
proximation method. It was first assumed that the surface can 
be described by a second-order equation and At = k bcl BC. If 
the deviations of the estimated temperature values from the 
experimental ones did not exceed 10°C, the calculations 
ended. With larger deviations, the assumptions of an equa¬ 
tion of the third, fourth, etc. order were investigated. 

In processing the available experimental data on the 
liquidus temperatures and the viscosity logarithms under a 
preset temperature, their relationship with the composition 
was approximated by polynomials of the nth degree. 

To assess the true value of the parameters taking into ac¬ 
count the required reliability level, it is possible to determine 
the confidence interval depending on the quantity of avail¬ 
able experimental points and the degree of the polynomial: 


bj ±tj 




where tj = t(P,k ) is the Student criterion for the number 
of degrees of freedom k = n—j— 1 [10]; S is the sum of 
the squares of the deviations of the experimental points v. 

n 

from the estimated points pj (x-); S =2^ l [y i -b 0 p 0 (x t )- 

i = l 
n 

biPi(Xi ) -... - b v p v (x t )] p) ( x i )• 

i = 1 

In this case all measurements are assumed to be equally 
precise and independent with normally distributed errors. 

Thus, the confidence interval is determined by the num¬ 
ber of degrees of freedom, which, in turn, depends on the 
number of experimental points and the degree of the polyno¬ 


and with the number of degrees of freedom equal to 9 

£,•±0.611 — , 

\ H j 

i.e., the interval increases 2 x 1,74 = 3.48 times with the 
number of degrees of freedom decreasing by 5 (9 - 4). 

With the confidence probability 0.95 these coefficients, 
respectively, are equal to 1.388 (with the number of the de¬ 
grees of freedom equal to 4) and 0.754 (the number of de¬ 
grees of freedom equal to 9), i.e., a decrease in the number of 
the degrees of freedom by 5 results in 2 x 1.84 = 3.68 times 
growth of the confidence interval. 

With a sufficiently large number of the experimental 
points, the degree of the polynomial can be increased to 
achieve a higher accuracy, and in the case of an insufficient 
quantity of the experimental data, it makes sense to restrict 
oneself to a linear dependence, since a decreased number of 
degrees of freedom will increase the confidence interval, i.e., 
will decrease the precision and increase the approximation 
error. 

Fluctuations in the chemical and mineral composition of 
the material can influence other important process parame¬ 
ters as well, in particular, the melt viscosity. The viscosity of 
melts and glasses has an effect on such important process pa¬ 
rameters as the temperatures of melting, clarification, work¬ 
ing, and annealing of glass. Furthermore, fluctuations in the 
composition lead to fluctuations of various significant pro¬ 
perties of the finished products, in particular, their TCLE and 
acid and alkali resistance. A variation in the composition can 
affect the liquation process and thus influence the decorative 
properties of the material. 

The authors have summarized data on the effect of varia¬ 
tions in the initial material composition on the properties of 
the batch, melt, and finished glass and glass ceramic materi¬ 
als in the quartz - albite (orthoclase, microcline) - anorthite - 
diopside - enstatite - magnetite systems corresponding to 
natural basalt, diabase, gabbro, norite, pyroxenite, and other 
basic and hyperbasic rocks. These data can serve as the basis 
for the mathematical description of property-versus-compo- 
sition dependences required for the automation of the calcu- 
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lations of such systems and the development of software pro¬ 
grams for fast control of the process parameters in produc¬ 
tion of glass and other silicate materials with prescribed 
properties. 

The Institute of Chemistry has developed a program for 
calculating the process parameters in the mineral fiber pro¬ 
duction, has tested this program at the Severonikel Integrated 
Works, and has obtained positive results. 
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